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ABSTRACT 

We explore the imprints of deviations from Gaussian primordial density fluctuations on the 
skeleton of the large-scale matter distribution as mapped through cosmological weak lensing. 
The skeleton of the Large Scale Structure traces its filamentary structure and thus provides a 
complementary piece of information with respect to, e.g., the simple Probability Distribution 
Function or the peak number counts. We computed the skeleton length of simulated effective 
convergence maps covering ~ 35 sq. deg each, extracted from a suite of cosmological 
n-body runs with different levels of local primordial non-Gaussianity. The latter is expected 
to alter the structure formation process with respect to the fiducial Gaussian scenario, and 
thus to leave a signature on the cosmic web. We found that alterations of the initial conditions 
consistently modify both the cumulative and the differential skeleton length, although the 
effect is generically smaller than the cosmic variance and depends on the smoothing of the 
map prior to the skeleton computation. Nevertheless, the qualitative shape of these deviations 
is rather similar to their primordial counterparts, implying that skeleton statistics retain good 
memory of the initial conditions. We performed a statistical analysis in order to find out at 
what Confidence Level primordial non-Gaussianity could be constrained by the skeleton 
test on cosmic shear maps of the size we adopted. At 68.3% Confidence Level we found 
an error on the measured level of primordial non-Gaussianity of A/nl ~ 300, while at 90% 
Confidence Level it is of A/nl ~ 500. In both cases we ignored the effects of observational 
noise and degeneracy with other cosmological parameters. While these values by themselves 
are not competitive with the current constraints, weak lensing maps larger than those used 
here would have a smaller field-to-field variance, and thus would likely lead to tighter 
constraints. A rough estimate indicates A/nl ~ a few tens at 68.3% Confidence Level for an 
all-sky weak lensing survey. 

Key words: cosmology: theory - gravitational lensing: weak - cosmological parame- 
ters - large-scale structure of the Universe 



1 INTRODUCTION 



One of the most striking successes of the standard cosmological 
paradigm consists in explaining the structure formation process 
in the Universe. Amplification of quantum fluctuations during the 
early (inflationary) phase of expansion led to seed density perturba- 
tions in the dark matter Overdense regions subsequently grew up 
due to gravitational instability, and eventually collapsed, detaching 
from the overall expansion of the Universe and forming the bound 
dark matter halos that today host galaxies and galaxy clusters. The 
standard model of inflation ( IStarobinskii|[T97i lGutlj|l98ll : iLindd 
1 1982 ) predicts a distribution of primordial density fluctuations that 
is indistinguishable from a Gaussian. This, coupled with the re- 
sults of perturbation theory and cosmological numerical simula- 



tions, gives us a relatively complete picture for the formation of the 
Large Scale Structure (LSS henceforth). 

Current cosmological experiments however still leave room 
for alternative scenarios. Particularly, extensions of the most stan- 
dard model of inflation can produce substantial deviations from 
a Gaussian d istribution of primordial densit y and p otential fluc- 
tuations (see lBartolo et al.ll2004lCheij|201(]| : lDesiacaues & Selii^ 
l20ld for recent reviews). The amount and shape of these deviations 
depend critically on the kind of non-standard inflationary model 
that one has in mind. Primordial non-Gaussianity is expected to 
modify the low-redshift LSS in different ways, most notably by 
shifting the timing for the collapse of the largest dark matter halos 
and inducing a model-specific scale dependence on the large scale 
bias thereof. Turning the argument around, observations of the LSS 
are critically important in order to validate or discard inflationary 
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models, and thus get a better grasp on the physics of the early Uni- 
verse. 

As a matter of fact, the problem of constraining devi- 
ations from primordial Gaussianity by means different than 
the Cosmic Microwave Background (CMB) intrinsic tem- 
perature anisotropies has recently attracted much attention 
in the literature, with effo rts directed towards the abundance 
of n on-linear s tructures ([M atarrese, Verde, & Jimenez 2000|: 
Verde et al.. 20001; iMathis. Diego, & Silk 2004: Grossi et al. 2001. 
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Maggiore & Riottd 12^01), halo biasing (e.g., iDalal et al 
'McDonaldl \mm . iFedeli, Moscardini, & Matarresd 



galaxy bispe ctrum (see 



for 

Ijeong & KomatstJ 



instan ce, ISefusatti & KomatstJ l2007l: 

l2009l) . mass density di stribution jGrossi et al.l |2008|) , cos- 
mic shear jpedeli & M oscardini! 1201 Ot IPace et al.l l201ll : 
iMaturi, Fe deli. & Moscardini 2011), integrated Sachs-Wolfe 
effect jAfshordi & ToUev 2008; Carbone, Verde, & Matarrese 
I2OO8I). Lv tt flux from low-density interg alactic medium 



(IVieletal.1 FoOa) 



21 -cm fluctuations dCooravl I2OO6I ; 

IPillepich, Porciani. & Matarrese! l2007h and reionization 
dCrociani et al."2009l). 

The gravitational deflection of light is one of the most pow- 
erful tools to probe the LSS of the Universe. Gravitational lensing 
reacts only to the total gravitational potential integrated along the 
line of sight, thus it is insensitive to the intrinsic nature and dynam- 
ical status of the matter providing the deflection. In cosmological 
context it can be used to probe the large-scale matter distribution 
without relying on strong assumptions, as instead is the case for bi- 
ased tracers such as galaxies and galaxy clusters. In the present pa- 
per we were interested in evaluating the effects of primordial non- 
Gaussianity on the topology of the (projected) matter distribution 
measured through cosmological weak lensing (cosmic shear). 

Topological statistics such a s the Minkowski F unction- 
als of the LSS ("Matsubara '2003'; 'Hikage et al.' '2008") and of 
the CMB maps ( Hikage, Komatsu, & Matsubara 2006) were al- 
ready recognized as useful tools fo r constraining the cos- 



mological initial conditions (see also iMatsubara & Jainl 1200 ll; 
VitelU. Jain. & Kamienll20()9l and the early results of lMessina et al.l 



I99(]| and lMoscardini et alj|l99 l'). In this work we dealt with one 



particular morphological indicator of cosmic shear maps, namely 
the skeleton of the effective convergence. The skeleton is a set 
of connected lines that trace the filamentary structure of a given 
field, and the rigorous theory for the skeleton of a two and three- 
dimensio nal random field can be foun d in Pogosvan et al. (2009D 
(see also iNovikov. Colombi. & Dord I2OO6; Pich on et a l. 201^. 
Models with a positive level of primordial non-Gaussianity are ex- 
pected to provide an anticipated structure formation, with the con- 
sequence that at early times the filamentary structure is presumed to 
be more pronounced with respect to models with a negative skew- 
ness. At late times however, the matter is also expected to flow 
more efficiently toward the knots of the cosmic web, thus leaving 
the filaments relatively emptier in the former models compared to 
the latter. Understanding the net effect on the filamentary structure 
of cosmic shear maps is an interesting issue by itself. 

The remainder of this work is structured as follows. In Section 
|2]we describe the basics of cosmologies with non-Gaussian initial 
conditions and the n-body/ray-tracing simulations that have been 
used for the subsequent analysis. In Section |3] we summarize how 
the skeleton of a two-dimensional field can be computed, and the 
relative local approximation that we used. In Section|4]we describe 
our results on the skeleton length in non-Gaussian cosmologies, 
and in Section[5]we summarize our conclusions. 



2 PRIMORDIAL NON-GAUSSIANITY 
2.1 Cosmological models 

A particularly convenient (although not unique) way to describe 
generic deviations from Gaussian initial conditions consists in writ- 
ing the gauge-invariant Bardeen's potential O as the sum of a 
Gaus s ian random field and a quadratic correct ion (Salopek & Bondl 
I99(]| ; lGanguietalJll994 I Verde et al.ll2OO0l ; iKomatsu & Spergel 



20011) . according to 



« = <I>G + /nl * (f I - <*G>) • (1) 

The parameter /nl in Eq. ^ determines the amplitude of non- 
Gaussianity, and it is in general dependent on the scale. The 
symbol * denotes convolution between functions, and reduces to 
standard multiplication upon constancy of /nl. In the following 
we adopted the large-scale structure convention (as opp osed to 
the CMB convention, see Afshord i & Tplley 2008; Car bone et al.l 
I2OO8I ; IPillepich, Porciani, & Hahnii2009l and lCrossi eTal]|2009h for 
defining the fundamental parameter /nl- According to this, the pri- 
mordial value of O has to be linearly extrapolated at z = 0, and 
as a consequence the constraints given on /nl by the CMB have 
to be raised b y ~ 30 per cent to comply with this pa per's conven- 
tion (see also iFedeli. Moscardini. & Matarresdl2009l for a concise 
explanation). 

In the case in which /nl the potential O is a random 
field with a non-Gaussian probability distribution. Therefore, the 
field itself cannot be described by the power spectrum Pt,(k) = 
Bld'^* alone, rather higher order moments are needed. In many cir- 
cumstances the dominant higher order moment is the bispectrum 
Bf,(ki , k2, Aa), that is the Fourier transform of the three-point cor- 
relation function (<l)(Ai)0(^2)'l'(^3)). It can be implicitly defined 
as 

<<D(*i)<D(*2)(l)(A3)> = (2nfSo(ki+k2 + k,)Bt,(ki,k2,k3). (2) 

As mentioned above understanding the shape of non- 
Gaussianity is of fundamental importance in order to pinpoint the 
physics of the early Universe and the evolution of the inflaton field 
in particular. The standard single-field inflationary scenario gen- 
erates negligibly small deviations from Gaussianity. These devia- 
tions are said to be of the local shape, and the related bispectrum of 
the Bardeen's potential is maximized for squeezed configurations, 
where one of the three wavevectors has much smaller magnitude 
than the other two. In this case the parameter /nl must be a con- 
stant, and it is expected to be of t he same order of the slow-roll pa- 
rameters jpalk, Rangaraian. & Srednickii l993) . that are very close 
to zero. 

However non-Gaussianities of the local shape can also be gen- 
erated in the case in which an additional light scalar field, differ- 
ent from the inflaton, contributes to the observed curvature pertur- 



bations teabich, Creminelli, & Zaldarriagal |2004|). This happens 
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for instance, in curvaton models ( Sasaki. Valiviita, & Wandsl2006l: 



'Assadullahi. Valiviita. & Wands' 2007 |) or in m ulti-fields models 
(Bartolo, Matarrese, & Riotto 2002; B emardeau & Uzan 2002). In 
this case the parameter /nl is allowed to be substantially different 
from zero, and the bispectrum of the primordial potential assumes 
the simple form 

Bt,{k, ,k2,ki) = 2fyi,B^ [l<:rX-^ + k"cX-^ + l<!r%'*] ■ (3) 

As naively expected, the mass function of dark matter halos 
gets modified in non-Gaussian cosmologies. Several prescriptions 
have bee n developed over the years in order t o describe this modi- 
fication jMatarrese. Verde. & Jimenezllioool ILo Verde et"al1l2008l; 
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iD'Amico et alj 1201 and all of them are in broa d agreement 
with numerical simulations of structure formation dOrossi et alJ 
l2009l : lwagner. Verde. & Boubekeuj2 01Cf). Only quite recently was 
it recognized that primordial non-Gaussianity introdu ces a scale 
dependence on the linear bias of dark matter halos jPalal et all 
l2008l : iMatarrese & Verdell2008l) . Such a dependence is regulated 
by suitable integrals of the potential bispectrum, and is one of 
the most promising LSS p robes of primordial non-Gau s sianit 



jAfshordi & ToUev l2008l: ICarbone, Verde. & Matarresel 



IVerde & Matarresei2009l;ICarbone. Mena. & Verdell2010l) . 



is siamty 

] llooi 



2.2 Numerical simulations 

The ;j-body simulations with non-Gaussian initial condit i ons tha t 
we employe d are origin ally described in iGrossi et al.l (l2007h : 
iHikage et al] (2OO8); Gr ossi et all ( I2OO8I) . We thus refer to those 
papers for detailed information, limiting ourselves to describe 
here only the main features. Dark matter-only structure formation 
was simulated by using the publicly available GADGET 2 code 
jSprin gel' 2005). assuming a flat ACDM cosmological model with 
density parameters for matter, baryons and cosmological constant 
equal to fi„,o = 0.30, fib.o = 0.04, and = 0.70, respectively. 
The Hubble constant was set to Hf, = A 100 km s"' Mpc"' with 
h = 0.7, while the matter power spectrum was normalized accord- 
ing to the constraint erg = 0.9. The comoving side length of the sim- 
ulation box measures 500/i"' Mpc, and the number of particles used 
is 800', implying a mass resolution of 2x IO'^/j^'Mq. In addition to 
the Gaussian case /nl = 0, six additional non-Gaussian scenarios 
were simulated, with the local shape for the primordial bispectrum 
described in the previous subsection and /nl = ±100, +500, and 
±1000. It should be noted that current constraints on /nl arising 
from the CMB and the LSS for the local case are A/nl < 100, thus 
only the non-Gaussian simulations with /nl = ±100 can be consid- 
ered realistic. Nevertheless, we performed our analysis for all the 
/nl values at our disposal both for illustrative purposes and because 
this allows us to better interpolate the results for other values. 

The ray-tracing analysis of the simulation cubes and the sub- 
sequent productio n of light cones and cosmic shear maps is de- 
scribed in detail in lPace et al. l|). Again, we merely summarize 
the procedure here. As a first step, the particles in the simulation 
boxes at each redshift snapshot have been randomly shifted and ro- 
tated in order to avoid repetitions of the same structure at different 
redshifts. Then, the particles themselves have then been projected 
along the line of sight on a plane lying in the middle of the box 
in order to compute the surface mass density on a regular grid. By 
the use of the two-dimensional Poisson equation this surface mass 
density has been subsequently converted into the lensing potential 
for the considered lens plane, from which it is possible to piecewise 
reconstruct the trajectory of each light ray shot from the observer 
toward the sources. Here all sources have been placed at the same 
redshift, Zs = 4. Although this might be not entirely realistic, we 
are here just interested in the relative differences between Gaussian 
and non-Gaussian cosmologies. The relation between the true and 
apparent positions of sources allows to compute a variety of lensing 
observables, including shear, flexion, and effective convergence. 

In this work we focused on the effective convergence only for 
two reasons: i) it is more directly related to the projected matter dis- 
tribution than, e.g., the shear, and thus it should have a more direct 
link to the initial conditions; ii) as we verified, the corresponding 
skeleton tends to be more regular than the skeleton of other weak 
lensing observables. The (two-dimensional) effective convergence 



can indeed be related to the three-dimensional matter distribution 

as 
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dw 



fK(w)fK(w^ - W) S [fK(w)e, w] 



(4) 



fxiws) a(w) 

where c is the speed of light, w = w(z) is the comoving distance 
out to redshift z, a(w) is the corresponding scale factor, d(x) is 
the (dark) matter density contrast at a given spatial position x, and 
fxiws - w) is the comoving angular diameter distance between the 
redshift z corresponding to w(z) and the source redshift, depending 
on the spatial curvature K (which is zero in our case). 

From each of the seven n-body cosmological simulations de- 
scribed above we extracted 60 effective convergence maps, each 
one obtained with a different randomization of the particle posi- 
tions in the simulation boxes. They can thus be regarded as 60 in- 
dependent realizations of the same light cone, allowing a mean to 
estimate the effect of cosmic variance. Each map has a resolution 
of 2048^ gridpoints, and has an opening angle of ~ 5.9 deg, corre- 
sponding to a covered area of ~ 35 sq. deg. 



3 THE SKELETON OF A TWO-DIMENSIONAL FIELD 

The skeleton of a two-dimensional field p is visually in- 
tended as the set of lines that trace the filamentary structure 
of the field itself. This naive de finition can be made rigorous 
jNovikov, Colombi. & Dorel2006l) by saying that the skeleton lines 
are parametrized by those trajectories x(t) = [xi(t), X2(t)] on the 
plane that obey the equation of motion 



— = Vp, 
dt ^ 



(5) 



starting at a saddle point and with initial velocity parallel to the 
direction of minimum curvature. This definition is non-local, and 
therefore difficult to be implemented practically. Fortunately, an 
acceptable local approximation of the s keleton can be devised 
jPogosvan et al.ll2009l ; |Pichon et aUbOld) . This local skeleton can 
be defined as a subset of the critical lines of the field, i.e. lines 
where the gradient of the field p is parallel to one of the eigenvec- 
tors of the Hessian matrix. If A is an arbitrary eigenvalue of the 
Hessian Hn, then this condition translates into 



■KpVp = AVp. 



(6) 



By adopting the index notation, Eq. ^ can be rewritten in a 
more useful fom as 



PijPj = Api, 



(7) 



where p; = dp/dxj, Pij = d'^p/dxjdxj, and summation over repeated 
indices is implicit. The two conditions in Eq. (O can be rewritten 
in a more insightful way, as 



Pii +P12 — = A, 
Pi 



P22 +P12 



P2 



(8) 



(9) 



By subtracting Eq. l|9j from Eq. (|8) we obtain after some manipu- 
lation that 

S =pipi{pii-P2i)+Pi2{pl-pi) = 0. (10) 

Note that Eq. dlOl l is valid whenever the gradient of the field p 
is parallel to one eigenvector of the Hessian, irrespective of which 
eigenvalue we are refeiTing to. The relevant eigenvalue is specified 
only in the pair of Equations ([8} and ([9}. The same Eq. JlOt can 
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Figure 1. The effective convergence map of a randomly chosen Hght cone 
reahzation for the Gaussian simulation (/nl = 0). smoothed with a Gaus- 
sian kernel whose rms is f = 0.02 times the side length of the map. The 
oveiimposed white lines represent the local skeleton of the smoothed field, 
as defined in the text. The side length measures ~ 5.9 deg. 



also be obtained along a different route, by simply imposing the 
condition that the gradient of the field is parallel to an arbitrary 
eigenvalue of the Hessian, i.e.. 



(11) 



where Ejj is the fully antisymmetric two-dimensional Levi-Civita 
tensor. 

The local skeleton can now be defined as the subset of the crit- 
ical lines for which A is the largest eigenvalues, i.e. the eigenvalue 
related to the direction of minimal curvature at a given point, and 
the smallest eigenvalue is negative. If we label the two eigenval- 
ues of the Hessian matrix as X\ and X2, and assume without loss 
of generality that Ai ^ A2, then the local skeleton is defined by 
Eqs. J8]l and ^ by setting A = Ai, with the supplemental condition 
that /li -I- /I2 ^ 0. As can be easily seen, since Ai ^ A2, the latter 
condition implies that A2 ^ 0, irrespective of Ai. 

In order to ensure the differentiability of the field under con- 
sideration up to second order, it is custumary to smooth the field 
itself prior to the calculation of the local skeleton. It is immedi- 
ate to understand that the length of the skeleton will depend on 
the smoothing length, with smoother fields having less features, 
and thus a shorter skeleton, with respect to less smooth ones. In 
this work we adopted a Gaussian smoothing kernel, and performed 
calculations for two different smoothing scales: C = 0.02 and 
£ = 0.002 times the side length of the maps. These values cor- 
respond to { ~ 1 arcmin and f ~ 43 arcsec, respectively. In the 
remainder of this work we will show alternatively the results ob- 
tained with the first or with the second smoothing scale, depending 
on which one fits best the purpose of illustration. 

As a demonstration, in Figure[T]we show the local skeleton of 
the effective convergence for one of the 60 light cone realizations 
in the reference Gaussian cosmology. The map has been smoothed 
with a Gaussian kernel whose rms i&l = 0.02 times the side length 
of the box. We observe that, since the smoothing has been per- 
formed with a Fast Fourier Transform (FFT) algorithm, the result- 
ing smoothed map has periodic boundaries. We also note two addi- 
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Figure 2. The skeleton length for excursion maps constructed from the ef- 
fective convergence map shown in Figure [T](/iml = 0), as a function of the 
threshold value in units of the field rrns. The skeleton length is expressed in 
units of the side length of the map (~ 5.9 deg), and five different sizes of 
the Gaussian smoothing kernel are adopted, as labeled in the plot. 



tional facts: i) the local skeleton follows quite nicely the ridges of 
the effective convergence map, highlighting the filamentary struc- 
ture thereof; ii) the local skeleton itself is quite irregular, and often 
broken in many pieces. This is a consequence of the local approx- 
imation, since the true skeleton w ould be made of continuos lines 
( iNovikov. Colombi. & Doril2006l) . 

The main statistics we will be interested in for the remainder 
of this paper is the length of the local skeleton related to excursion 
maps. Excursion maps are subsets of the original map given by 
considering only those points where the field value exceed some 
threshold value, p ^ p,. For subsequent use, we defined a new 
field 77 that is given by the original field divided by its rms, that 
is 7/ = pier, where a = ^j{p-) - {pY- Under this rescaling the 
threshold field value maps into a threshold rj value, = pja. As 
in Figure[T]and in the remainder of this work, p = k^, where is 
the effective lensing convergence. We stress that, since the effective 
convergence is an integral along the line of sight of the density 
contrast, it may also take negative values. 

Measuring the length of the skeleton i s an issue that de- 
serves some further discussion. In other work s ( Eriksen et alj2004l : 
iNovikov. Colombi, & Doril2006l : IHou et al.ll20I(]|) this length has 
been measured by using linear or cubic interpolation in order to 
find where a skeleton line effectively crosses the borders of a pixel 
that is known to host a critical curve. Then, it was assumed that the 
curve itself is a line segment between the crossing points on two 
sides of the pixel. Of course, this reasoning does not work when 
multiple critical lines converge to the same point (critical point) 
within one pixel or when the critical curves are interrupted, an in- 
stance that is often verified for the local approximation of the skele- 
ton at our resolution level (see Figure[T). Thus, since some level of 
approximation is unavoidable, we chose to adopt a simplified ap- 
proach, and assign to each pixel a length scale that equals its side 
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Figure 3. The distribution of tlie total skeleton length (in units of the map side length) for effective convergence maps smoothed with a Gaussian kernel of 
width C = 0.02 times the side length (~ 7 arcmin, left panel) and C = 0.002 times the side length (~ 43 arcsec, right panel). Distributions are computed over 
the 60 light cone realizations, and results for three different values of the non-Gaussian parameter /nl are shown, as labeled in the plots. Note that the abscissa 
in each panel is shifted according to the average value of the total skeleton length in the Gaussian case. 



length. The average length of a line segment randomly crossing a 
squared pixel is probably somehow different from the side length 
of the pixel. However, we are not interested in the absolute value of 
the skeleton length, rather on how this and its differential and cu- 
mulative counterparts are modified in presence of primordial non- 
Gaussianity. Thus, as long as we adopt the same convention for all 
maps in all simulations, our conclusions should hold. 

In Figure |2] we show, for the same effective convergence re- 
alization represented in Figure [T] the length of the local skeleton 
for excursion maps as a function of the threshold z/t used to define 
the excursion map itself. We name this quantity as the cumulative 
skeleton length -£(?7t). while the differential skeleton length is its 
derivative normalized by the length of the skeleton of the original 
map, D (r]t) = - [l/JK-oa)] d£,(r]t)/dr]f Another quantity we shall 
be interested with is the normalized cumulative skeleton length, 
namely C(77t) = £.{rit)l £.{-oa). In Figure |2] we report curves for 
several different values of the Gaussian smoothing scale around 
I = 0.02. As one could naively expect, increasing the smoothing 
scale decreases the length of the skeleton, since the more the field 
is smooth, the less features it displays. 



4 RESULTS 

4.1 Total skeleton length 

We now turn attention on how the skeleton is modified by the intro- 
duction of some level of primordial non-Gaussianity. The first thing 
that we have investigated is how the total length of the local skele- 
ton (without excursion) does change. It is expected that primordial 
non-Gaussianity with a positive skewness of the matter density field 
would produce an early enhancement of the filamentary structure in 
the large scale matter distribution, that would be later compensated 
by the more efficient growth of dark matter clumps. On the other 
hand, a negative skewness would produce the opposite effect. 



In the left panel of Figure |3] we show the distribution of the 
total length of the local skeleton -£(-oo) obtained by the 60 effec- 
tive convergence maps we studied. In order to produce this panel 
all maps have been smoothed with a Gaussian kernel whose rms is 
f = 0.02 times the side length of the map (corresponding to ~ 7 ar- 
cmin). Besides result for the fiducial Gaussian cosmology, we also 
report those for the most extreme non-Gaussian simulations, with 
/nl = +1000. Despite having adopted the largest /nl values at our 
disposal, it is apparent that very little difference is present between 
the three distributions, which overlap with each other almost com- 
pletely. It is true that the average value of the distribution referring 
to /nl = +1000 is larger than that relative to /nl = -1000, im- 
plying an effective net enhancement of the filamentary structure, 
however the difference between the two is only ~ 0.25 (times the 
side length of the map), while the rms of the individual distributions 
is ~ 0.6. 

In the right panel of Figure |3] we show the same results, but 
this time assuming a smaller smoothing length of ^ = 0.002 times 
the side of the map (corresponding to ~ 43 arcsec), in order to 
check whether a larger difference in the total skeleton lengths could 
arise upon inclusion of smaller scale features of the effective con- 
vergence maps. As it turns out, this is indeed the case: the three 
distributions show now a substantial separation in the same direc- 
tion noted above. In this case the average values of the two distri- 
butions referring to the cases /nl = ±1000 differ by ~ 1.6 times 
the side of the map, while the rms is ~ 0.8 in the same units. This 
result is expected, and can be explained as follows. Primordial non- 
Gaussianity tend to affect more strongly the most extreme features 
of the LSS, for instance the high-density peaks and ridges. The ef- 
fective convergence we are dealing with is actually a projection of 
the LSS, and the larger is the smoothing radius, the more these 
extreme features are washed out from the maps, thus effectively 
reducing the leverage of the skeleton to tell the models apart. A 
simple Kolmogorov-Smimov test implies the two distributions for 
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Figure 4. The differential skeleton length distiibution (black solid line) and the PDF (red dashed line) of the effective convergence for the Gaussian cosmology. 
The lines show the average values computed over the 60 light cone realizations as labeled, while the shaded regions represent the dispersions around the means 
of the distributions for each j/f Each effective convergence map has been smoothed with a Gaussian kernel of radius t = 0.02 (left panel) and ( = 0.002 times 
its side length (right panel). 



/nl = -1000 and /nl = +1000 being the same at more than 20% 
Confidence Level for £ = 0.02 and only at less than 0.5% Con- 
fidence Level for £ = 0.002, further confirming the result of the 
visual inspection. 

The fact that the dispersions of the three distributions shown 
in Figure [3] is larger for £ = 0.002 than for £ = 0.02 does not 
mean that the total skeleton length is more affected by the cosmic 
variance in the former case, because is much larger 

there than in the latter. As we will detail in the next Subsection, the 
normalized cumulative and differential skeleton lengths, that are 
correctly normalized, actually show the opposite behavior. 



4.2 Differential and cumulative skeleton lengths 

The next step is to deal with the differential length of the lo- 
cal skeleton D(t]i) and its cu mulative counterpart C(?7t). We star t 
with the former. As noted in iNovikov. ColombiT & Dorel tOOdt) 
the differential skeleton length should behave quite similarly to the 
Probability Distribution Function (PDF henceforth) of the field, al- 
though several subtle differences exist that make it worth study- 
ing JHou et alteOlOl ). One such difference is made by the fact that 
we are considering the skeleton of the effective convergence, but 
not the anti-skeleton (that would follow the valleys of the field in- 
stead of the rid ges, and can be viewed as the skeleton of -k^, see 
IPogosvan et al . 2009). As a consequence the differential length dis- 
tribution is expected to be more skewed toward high effective con- 
vergence values as compared to the simple PDF. This is an advan- 
tageous feature for our purposes, since primordial non-Gaussianity 
has indeed more effect on the more extreme details (e.g., ridges) 
of the LSS. Although not included here, the anti-skeleton by it- 
self could share equally advantageous features, since primordial 
non-Gaussianity is expected to affect the most extreme voids as 
well. Even if the anti-skeleton was included in the analysis how- 
ever, the critical lines would still sample the field in a non-random 



way, at least if we are not in the peculiar, completely incoherent 
case. Thus, difference s between the PDF and the differen tial length 
are expected anyway (iNovikov. Colombi, & Doril2006h . 

In Figure|4]we report the differential length distribution of the 
skeleton D(r]t), comparing it with the PDF of effective convergence 
maps in the reference Gaussian cosmological model. As expected, 
the differential length of the local skeleton is more shifted to high 
values of the effective convergence with respect to the PDF. This 
is in agreement with the naive expectations that we outlined above. 
Particularly, the peak of the differential length is at ~ 0.5cr, while 
the peak of the PDF is only at 0.5it. In the same Figure we re- 
port results for the two different smoothing lengths adopted in this 
work, £ = 0.02 and £ = 0.002 times the side length of the maps. 
The behaviors of the distributions in the two cases are relatively 
similar. Specifically, the positions of the peaks are approximately 
unchanged. The only significant difference is given by the disper- 
sion around the mean for a fixed r]t, which is substantially larger 
in the case £ = 0.02 than in the case £ = 0.002. This behavior is 
expected because for smoother fields the typical length scale of fea- 
tures tend to be larger, thus they are presumed to be more subject 
to field-to-field variance. 

An alternative way to display the difference between the PDF 
of the effective convergence field and its differential skeleton length 
can be seen in Figure|5] where we show the distribution of the ratio 
of D(77i) to PDF(77i) in the fiducial Gaussian cosmological model. 
As a consequence of the tendency of the local skeleton to trace 
the ridges of the map, this ratio increases steadily with increas- 
ing excursion threshold The trend of the rat i o is q ualitatively 
the same that has been shown in lPogosvan et alj i2009l) for a two- 
dimensional Gaussian random field, although the field we are con- 
sidering is not Gaussian, even if the initial conditions are, due to 
the non-linear clustering of dark matter. 

The normalized cumulative skeleton length distribution for the 
fiducial Gaussian cosmology is shown in Figure|6]for the two Gaus- 
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Figure 5. The distribution of the ratio of the differential skeleton length to 
the PDF of the effective convergence in the Gaussian simulation. For each 
value of rif the average value over the 60 light cone realizations is shown by 
the black solid line for C = 0.02 and by the red dashed line for I = 0.002, as 
labeled in the plot. The shaded regions show instead the dispersions around 
the mean values. 




-4 -2 2 4 

Figure 6. The normalized cumulative skeleton length distribution of the 
effective convergence for the Gaussian cosmology. The lines show the av- 
erage value over the 60 light cone realizations, while the shaded regions 
represent the dispersion around the mean of the distribution for each 
Each effective convergence map has been smoothed with a Gaussian kernel 
of radius t = 0.02 and C = 0.002 times its side length, as labeled in the plot. 



sian smoothing scales that we adopted in this work. The behavior 
is the one expected by looking at Figure |2] with only little differ- 
ence between the two smoothing radii adopted. Also in this case the 
only significant difference is given by the dispersion around the av- 
erage value, that is substantially larger when C = 0.02. As noted by 
iHouetalJjioTol) , the differential and normalized cumulative skele- 
ton lengths are equivalent quantities, and should return comparable 
information about the field. 

In order to visually compare the results for the fiducial Gaus- 
sian and non-Gaussian cosmologies it is customary to consider 
the distribution of the difference between the differential skeleton 
length of light cone realizations in the non-Gaussian models to the 
average value of the differential skeleton length in the Gaussian 
model. In other words, for each non-Gaussian cosmology the dis- 
tribution of the quantity D(?7t) - (X^'^^rji')^ is considered. A sim- 
ilar approach is followed for the normalized cumulative skeleton 
length, with the trivial replacement of 2)(?/t) with C(?7t). This is done 
because the differences between cosmologies that we would like to 
highlight are usually very small, such that they would be barely be 
visible in plots such as those shown in Figure|4] 

In Figure |7] we report the distribution of the difference be- 
tween the differential skeleton length for the non-Gaussian cos- 
mologies with /nl = ±1000 and the average value of the same 
quantity in the fiducial Gaussian model. We chose the two most 
extreme non-Gaussian models in order to better highlight the ef- 
fects. Also, we refer to the smoothing length ( = 0.002 times the 
map side length. As can be seen, the differential skeleton lengths 
in the two non-Gaussian simulations differ on average from their 
Gaussian counterpart of at most ~ 0.01, while the scatter around 
the mean in the Gaussian case can reach ~ 0.02. In other words, 
the effect of even a large level of primordial non-Gaussianity on 



the differential skeleton length would be substantially smaller than 
the field-to-field variance of the fiducial model. 

A similar conclusion can be reached by looking at Figure |8] 
where we show the same as in Figure |7] but for the cumulative 
length distribution. Again the effect of primordial non-Gaussianity 
tends to be significantly smaller than the cosmic variance estimated 
for the fiducial Gaussian cosmology. In both cases, we decided not 
to show the results for the larger smoothing length C = 0.02. These 
results show the same generic features than those for £ = 0.002, 
however the cosmic variance is much more important (see the 
discussion above). Hence we decided to stick just with the most 
promising case, to which we limit ou r subsequent analysis. 

In Figure 1 of IHou et alj ilOld) the authors show the equiv- 
alent of our Figures |7] and [8] for a perfectly Gaussian two- 
dimensional random field that has been perturbed by a non- 
Gaussian contribution with local shape. There are three main dif- 
ferences between their plots and ours: ;) the curves with positive 
and negative /nl are inverted, due to a sign difference in the def- 
inition of the differential and cumulative skeleton lengths (irrel- 
evant for our purposes); ii) our effective convergence maps can- 
not be a non-Gaussian random field with local shape, due to the 
non-li near growth of structures and the projection procedure; Hi) 
IHou e t al. (2010) do not make any distinction between the skeleton 
and the anti-skeleton. Despite these differences, it is remarkable 
that their plots look rather similar to ours. Specifically, the succes- 
sion of minima and maxima of the non-Gaussian curves and the 
generic shape of the dispersion around the mean for the Gaussian 
case are well cap tured. The main difference between the plots of 
IHou etaljj20ldh and ours resides in the fact that ours are not sym- 
metric around rjt = 0, and are more skewed toward higher values of 
the excursion threshold. This is a consequence of the combination 
of points ii) and Hi) above, in that the absence of the anti-skeleton 
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Figure 7. The difference between the differential skeleton length in the two 
non-Gaussian cosmologies with /nl = ±1000 and the average value of the 
same quantity in the reference Gaussian cosmology. The black solid line 
represents the average value of the distribution for the Gaussian cosmol- 
ogy (that is correctly zero), and the shaded region indicates the dispersion 
around the average. The other lines refer to the averages for the two non- 
Gaussian cosmologies, as labeled in the plot. Maps are smoothed with a 
Gaussian of radius t = 0.002 times their side length. 



Figure 8. The difference between the normalized cumulative skeleton 
length in the two non-Gaussian cosmologies with /nl = ±1000 and the 
average value of the same quantity in the reference Gaussian cosmology. 
The black solid line represents the average value of the distribution for the 
Gaussian cosmology (which is correctly zero), and the shaded region indi- 
cates the dispersion around the mean. The other lines refer to the averages 
of the non-Gaussian cosmologies, as labeled in the plot. 



in our analysis gives more importance to the ridges of the effec- 
tive convergence, while the non-linear clustering of dark-matter 
enhances the ridges themselves. Despite these differences, we can 
conclude that the skeleton of the effective convergence keeps a 
good memory of the statistics of initial density fluctuations, and for 
this reason it is a potentially powerful mean for probing primordial 
non-Gaussianity. 

4.3 Statistical analysis 

As a final step of our analysis, we performed a statistical investi- 
gation in order to understand what is the level of local primordial 
non-Gaussianity that could be detected by measuring the skeleton 
length of cosmic shear maps similar to those simulated here. The 
analysis is the same for both the differential length and the nor- 
malized cumulative length with just the replacement of 2)(?/t) with 
C{rit). Therefore, in the following we describe the formalism in the 
former case only. 

As a first step, we recall that the differential and normal- 
ized cumulative skeleton lengths have been computed by adopting 
y = 32 different threshold values rju, evenly spaced between -4 
and +A (we remind the reader that the threshold 77, is in units of 
the rms of k^). In the following discussion we adopt the notation 
2),(/nl) = Wnl,??!.,), and = 2),(0). Given this, we first of 
all computed the covariance between different thresholds for the 
Gaussian model, as 

T,, = (pf'Df) - ipf') (£)f ) , (12) 

where the average values are taken over the 60 light cone realiza- 
tions. 



Then, we assumed that the differential skeleton length of the 
fiducial Gaussian model was the measured one, with each light- 
cone realization representing an independent measurement, while 
the average value of the skeleton length for the various non- 
Gaussian cosmologies was the prediction of the assumed model. 
We thus defined a^-(/NL) function according to 

X\M = t [(n(/NL)) - 0;°'] [(2),(/nl)) - of] . (13) 

'7=1 

In principle, for a fixed value of /nl, one could compute ;^^^(/nl) for 
each of the 60 light-cone realizations, and then consider the average 
value. However, we verified that 60 realizations are not sufficient, 
in the sense that the average value of^^(O) does not converge to v 
as it is expected to. In order to compensate this, we adopted the fol- 
lowing procedure. For each pair (i, /) in the sum of Eq. l ll3t we uti- 
lized the values of and Of of a randomly chosen light-cone. 
We then repeated this procedure 10'' times (obviously allowing for 
repetitions) and considered the resulting average value of x^(/nl). 
Now, we correctly get;^^^(0) v = 32. 

We are therefore in the position to compute a Ax^ifNh) func- 
tion as Ax^ifhth) = ^"(/nl) - y, that vanishes for /nl = 0, and 
can thus be used to set confidence intervals on the parameter /nl- 
In Figure |9] we show the A;^'^(/nl) values computed for each of 
the simulations in our sample, and considering both the differen- 
tial and the normalized cumulative skeleton lengths. In order to 
set Confidence Levels (CL henceforth), we fitted the results with 
two parabolic curves. The fact that the A_;^'^(/nl) curves are well 
fitted by parabolae implies the interesting fact that both the differ- 
ential and the cumulative skeleton lengths are linear functions of 
/nl. Moreover, the curves resulting from the differential and the 
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Figure 9. The A;^'^(/nl) function for the differential skeleton length (red 
line and circles) and the nonnalized cumulative skeleton length (green line 
and triangles), assuming a Gaussian smoothing of effective convergence 
maps with radius f = 0.002 times their side lengths. The points represent 
the actual A^^ calculations, while the curves are the corresponding best 
parabolic fits. The hoiizontal dotted lines represent the values corre- 
sponding to various Confidence Levels, as labeled. 



cumulative skeleton lengths are very similar, coherently with the 
fact that the two estimators are equivalent and should return con- 
sistent results. Using standard A^^ variations corresponding to one 
parameter, it turns out that by computing the local skeleton of ef- 
fective convergence maps covering ~ 35 sq. deg would measure 
the level of primordial non-Gaussianity with an error A/nl ~ 300 
at 68.3% CL, and A/nl ~ 500 at 90% CL. As explained above, 
constraints obtained by maps with a smoothing length of C = 0.02 
would be considerably looser, due to the much larger impact of the 
cosmic variance. 



5 DISCUSSION AND CONCLUSIONS 

In this work we computed the local approximation of the skele- 
ton for simulated effective convergence maps. The maps are ex- 
tracted from cosmological simulations of structure formation with 
non-Gaussian initial conditions of the local shape. For each non- 
Gaussian cosmology we have 60 maps, that can be considered as 
independent realization of the same light cone, and each map cov- 
ers an area of ~ 35 sq. deg. The simulations span a wide range 
in the level of primordial non-Gaussianity /nl, having in particular 
/nl = 0, /nl = +100, /nl = +500, and /nl = +1000. In order 
to insure differentiability of the field we smoothed the maps with a 
Gaussian having two different length scales, namely ( = 0.02 and 
{ = 0.002 times the side length of the maps themselves. 

By considering excursion subsets of each map we evaluated 
the differential and cumulative length of the local skeleton, which 
carry equivalent cosmological information, yet alternative to more 
common statistics like the PDF or the peak number counts of the 



field. The main results that we derived in this paper can be summa- 
rized as follows. 

• The total skeleton length of effective convergence maps is 
larger in models with a positive /nl and smaller if /nl is negative. 
The separation between models with positive and negative /nl val- 
ues tends to vanish if the field is smoothed with a relatively large 
smoothing length. This is so because in this case the most promi- 
nent ridges of the field, that are presumably the most affected by 
primordial non-Gaussianity, tend to be washed out. 

• The result outlined in the previous point is in agreement with 
the fact that a positive-skewed primordial non-Gaussianity leads to 
an anticipated formation of structures, with the filamentary struc- 
ture being thus more pronounced. The subsequent flow of matter 
from the filaments to the knots of the cosmic web, which should 
also be enhanced in models with positive /ml, is not enough to 
counteract the first effect, at least at the redshifts probed by cos- 
mological weak tensing. 

• The distributions of the difference between the normalized cu- 
mulative and differential skeleton lengths and the average values of 
the same quantities in the fiducial Gaussian cosmology are qual- 
itatively rather similar to what one would expect for a Gaussian 
random field perturbed according to a local shape non-Gaussianity. 
This implies that these statistics retain good memory of the initial 
conditions, despite the subsequent non-linear growth of structures 
and the projection procedure. 

• Even if the most extreme non-Gaussian models are used, hav- 
ing /nl = ±1000, the average differences of the differential and 
normalized cumulative skeleton lengths to the corresponding av- 
erage values in the Gaussian case are always significantly smaller 
than the dispersion around the mean for Gaussian cosmologies. In 
other words, the signal we are trying to detect is always smaller 
than the noise given by the cosmic variance for the field size we 
have considered. This is more so for larger smoothing scales, since 
smoother maps are expectedly more affected by field-to-field vari- 
ations. 

• A simple statistical investigation based on intervals shows 
that the constraints expected by the skeleton analysis performed 
on cosmic shear maps similar to those used here are rather weak. 
Particularly, the expected error on the level of primordial non- 
Gaussianity is of A/nl ~ 300 at 68.3% CL and of only A/nl ~ 500 
at 90% CL. These numbers are substantially larger than the cuixent 
constraints given by the CMB and LSS. 

While the results concerning the constraints on /nl can be 
quite discouraging, it is worth stressing that they are valid only for 
the map size that has been considered here, ~ 35 sq. deg. Evidently, 
larger maps would be less affected by cosmic variance, and thus 
would lead to more stringent constraints on the level of primordial 
non-Gaussianity. For the same reason, it is safe to state that mea- 
suring the skeleton of the effective convergence in the COSMOS 
field (~ 2 sq. deg) or in the Deep Lens field (~ 4 sq. deg) would 
not lead to any signature of primordial non-Gaussianity. Although 
we are not in a position to properly quantify how much the covered 
area should be in order to get significant bounds o n /nl, we can per- 
form a roug h estimate along the following lines. In lSomerville et aU 
(2004) and lMoster et al.l ( 1201 ih the authors quantify how the cos- 
mic variance for dark matter density fluctuations scales with the 
volume of space that is probed. It turns out that cr^ oc V^^ with 
7 ~ 0.4 - 0.6, although a power law is probably not the best de- 
scription of the trend. For a fixed source redshift, the weak lensing- 
surveyed volume of space is proportional to the covered area, thus 
we can infer that the same scaling applies to the sky coverage as 
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well. Now, we are not looking at the dark matter density directly, 
rather at the skeleton of its projection. Still, given the similarities 
between, e.g., the differential skeleton length and the PDF of the 
effective convergence, we assume that the relative cosmic variance 
shares the same scaling with the survey area. This means that, for 
instance, the CFHTLS field, covering almost 5 times the area that 
we considered here (~ 170 sq. deg), would lead to A/nl ~ 200 at 
68.3% CL. In order to have an error on the estimated value of pri- 
mordial non-Gaussianity of A/nl ~ 100 at the same CL, we would 
need a sky coverage of ~ 300 sq. deg. Plan ned wide field cosmic 
shear surveys such as Euclid ( Laureiisl2009 ) will produce effective 
convergence maps of the entire extragalactic sky, thus providing a 
powerful test of the initial conditions by using the skeleton statis- 
tics, with an estimated A/nl ~ a few tens at 68.3% CL. 

Meanwhile, we would also like to outline a few limitations 
of our approach. We assumed perfect knowledge of the effective 
convergence in our maps, while real maps are affected by various 
sources of noise, both systematic and statistical, such as the intrin- 
sic ellipticity distribution of the sources used to measure the cos- 
mic shear. Furthermore, in our statistical analysis we were allowed 
to vaiy only the parameter /nl> while a more complete investiga- 
tion would consider also marginalizing over the other cosmologi- 
cal parameters (especially erg) that might affect the skeleton length. 
Performing this, and properly answering the above question about 
the minimum field size necessary to obtain competitive constraints 
on /nl would require a much larger and detailed set of simulations, 
that is not available at the moment. 

As this work was being finalized, another paper dealing 
with the topo l ogy of weak-lensing maps appeared on the archive, 
iMunshi et al.l ( l201ll ). That work is a nice complement to our, 
since it presents a perturbative approach in order to compute the 
Minkowski Functionals of the effective convergence. It further 
stresses the importance of investigating all pieces of information 
extractable from cosmological weak lensing (morphology, correla- 
tion function, peak counts, etc.), in order to maximize the return of 
present and future cosmic shear surveys. 
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